Characterization of complex quantum dynamics 
with a scalable NMR information processor 
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We present experimental results on the measurement of fidelity decay under contrasting system 
dynamics using a nuclear magnetic resonance quantum information processor. The measurements 
were performed by implementing a scalable circuit in the model of deterministic quantum compu- 
tation with only one quantum bit. The results show measurable differences between regular and 
complex behaviour and for complex dynamics are faithful to the expected theoretical decay rate. 
Moreover, we illustrate how the experimental method can be seen as an efficient way for either 
extracting coarse-grained information about the dynamics of a large system, or measuring the de- 
coherence rate from engineered environments. 
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Quantum information processors (QIP) promise effi- 
cient solutions to problems that seem intractable using 
classical devices. Tremendous progress has been achieved 
in understanding how much more powerful quantum me- 
chanics is at manipulating information Some first 
steps have been taken in experimentally demonstrating 
both algorithms Q and simulating physical systems 
in small prototype devices, and extensive research is 
pushing forward many avenues in the search for a fully 
scalable quantum computer Q. On the way to larger 
and scalable devices, benchmarking and characterizing 
will play a crucial role in understanding the performance 
of devices and pointing out directions where improvement 
is required. An obvious characterizing procedure is full 
state or process tomography [j| ; however, these methods 
will become unfeasible as we move to larger QIP. The 
number of experiments necessary for tomography grows 
exponentially with the size of the system and the exper- 
imental time required will quickly become prohibitive. 
We must turn instead, to characterizations which pro- 
vide coarse-grained information. These experiments will 
extract the information necessary for quantum control 
or error correction processes, yet will be achievable us- 
ing a reasonable amount of resources - i.e. in a scalable 
manner. 

One such measure of interest is the fidelity decay. This 
tool is useful in two contexts: characterizing the com- 
plexity of a system's dynamics; and, in studying the ef- 
fects of the environment dynamics on decoherence rates. 
Simulating quantum dynamics, whether of the system of 
interest or of a large environment, on a classical com- 
puter is inherently difficult; so, a quantum computer is 
a natural tool in the study of complex (chaotic) quan- 
tum dynamics [jj. Previous experimental work on the 
simulation of quantum chaos on a QIP relied on full 
state tomography, which is an inefficient, and hence un- 
scalable means of measuring signatures of the system's 



complexity. Here, we experimentally investigate complex 
dynamics through a fidelity decay measurement that re- 
quires only one pseudo-pure bit of quantum information, 
is fully scalable in the absence of noise, and whose run- 
ning time is limited only by the inherent complexity of 
the map we are characterizing. 

In classical mechanics, the distinction between regu- 
lar and chaotic behaviour is clear: the rate at which two 
states, initially close in phase space, are driven apart, 
is either exponential or sub-exponential. This measure 
has no relevance for isolated quantum systems, since 
their unitary evolution preserves the distance between 
two states. The field of quantum chaos looks at other sig- 
natures of chaos yj- Fidelity decay under perturbation, 
initially proposed by Peres |9(|, is an important analogue 
to the classical notion of instability, and has received 
much study 0, 0, 0, 0|. Rather than considering 
the overlap of two nearby states evolving under the same 
evolution, this method considers the overlap of the same 
initial state evolved under two nearby unitary maps. If 
the discrete time evolution is given by the map U , then 
we also consider the map U p = UP where P = exp (— iV) 
for some hermitian matrix V. Then, the fidelity decay 
after n steps is given by, 
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The analogy with classical mechanics suggests that 
chaotic systems might show an exponential fidelity de- 
cay, while regular systems would decay at a slower (e.g. 
polynomial) rate. However, numerical studies of small 
systems show a much more complicated situation, with 
different behaviours de p end ing on the form and strength 
of the perturbation [llLlia . Il3j. Yet, for a range of suffi- 
ciently strong perturbations, chaotic systems show a uni- 
versal exponential fidelity decay at an average rate which 
is given by the Fermi Golden Rule (FOR) tHHl- Thc 
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average decay is universal, in the sense that it is indepen- 
dent of the system dynamics and depends only on certain 
details of the perturbation. Specifically, the rate of the 
universal exponential decay depends only on the variance 
of the eigenvalues of the perturbation (i.e., the 2-norm of 
the perturbation, ||V||2) and continues until a saturation 
level O(jj), due to the finite system dimension N. Hence, 
direct measurement of the fidelity decay under different 
applied perturbations can provide important information 
about the complexity of the system dynamics. 

Given some dynamical system under study, the fidelity 
decay measured for any one particular initial state will 
fluctuate from the average - an effect particularly pro- 
nounced with small systems. Hence estimating the aver- 
age decay can require averagin g over many different ini- 
tial states. However, recently IPoulin et alJ ljj worked 
out a quantum circuit (see Fig. for directly measuring 
the average fidelity decay within the deterministic quan- 
tum computation with a sing lc bit (DQC1) model 0. 
The fidelity decay after n steps, averaged over a uniform 
(Haar) measure of the initial states \ip) (sec Eq. for an 
arbitrary system of dimension N takes the general form, 
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Therefore, experimental determination of the average 
fidelity decay for a particular U and P requires measuring 
the trace, for which there exists an efficient DQC1 circuit 
0. 
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FIG. 1: Ideal quantum circuit for measuring fidelity decay 
after n steps. The top qubit in a pure state can be considered 
either as a probe of the bottom system of K qubits which 
starts in the maximally mixed state or as a toy system be- 
ing decohered by the maximally mixed environment below. 
Making the final rotation about the x (y ) basis gives the real 
(imaginary) part of the trace in Eq. [5] Taken from Ref. [lEt 

Liquid state NMR offers a test-bed technology with 
a sufficient number of qubits and control for the first 
demonstration of a DQC1 algorithm. The DQC1 model, 
although more limited than full quantum computation, is 
interesting, since it questions entanglement as the basic 
concept behind the extra-power of quantum computers 
[l8| . Although the model does not allow for error cor- 
rection, the required initial state is not a full multi-qubit 
pseudo-pure state as in conventional NMR QIP [l^ i but 
a single pseudo-pure qubit of the form p t = i^ll+e|0){0|, 
with all other qubits starting in the maximally mixed 



state. This state is very close to the thermal state of liq- 
uid state NMR and can be efficiently prepared from the 
thermal distribution. 

The experiments were implemented on a Bruker 
Avance 700 MHz spectrometer using the molecule trans- 
crotonic acid (shown in Fig. The four carbons were 
used as the system qubits and the spin \ selected part 
of the methyl group was the probe. This setup allowed 
the use of broadband rcfocussing pulses to decouple the 
probe qubit from the system qubits during the applica- 
tion of the random map U. H± and H2 were decoupled 
by placing them in a pseudo-pure state at the beginning 
of the pulse sequence. Measurement was performed by 
quadrature detection of the free induction decay of the 
methyl group which gives both (a x ) and {u y ). Further 
details of the molecule and experimental technique can 
be found in Ref. E3 . 
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FIG. 2: The molecule trans-crotonic acid and its natural 
Hamiltonian - diagonal terms are chemical shifts and off- 
diagonal elements coupling terms in Hz. The darkly shaded 
nuclei are hydrogen and the lightly, carbon. The unlabeled 
nuclei are oxygen whose natural abundance have zero spin. 
Hence, we can ignore their coupling with the rest of the 
molecule. The three hydrogen nuclei of the methyl group 
are magnetically equivalent. They form a composite spin, 
from which we can select the spin | subspace, giving us an 
additional qubit. 

We investigated the difference in fidelity decay re- 
sponse for regular evolution - provided by the natural 
Hamiltonian of the molecule - versus complex evolution 
- achieved by using our ability to generate arbitrary 
dynamics (universal control) to implement a pseudo- 
random map 0. We used a decomposition for pseudo- 
random operators based on repeated applications of a 
two step process: (1) Individual random rotations are 
applied to each qubit; (2) Simultaneous two-body inter- 
actions between neighbours given by the unitary, U = 



exp 



|2jJ |. For a finite number of 

repetitions of these two steps, the resulting distribution 
of maps is biased with respect to the uniform measure; 
however, it converges exponentially to the uniform mea- 
sure . For this experiment the first step was replaced 
by correlated random rotations, induced by broadband 
pulses of random phase and power on either side of a ran- 
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domly chosen delay period providing random rotations 
about the z-axis and some two body interactions. In or- 
der to reach a suitable trade-off between decoherence ef- 
fects and the randomness of the map four iterations were 
implemented. Although there is not as many random pa- 
rameters in our experimental scheme as in the theoretical 
procedure, numerical simulations showed that, after four 
repetitions it implemented a sufficiently random map to 
simulate complex dynamics. The first group to imple- 
ment these pseudo-random operators measured their de- 
gree of randomness through complete state tomography 
[21|. Note that here, from the measurement of the fi- 
delity decay we can efficiently characterize the degree of 
randomness of the map. 

The controlled perturbation was implemented using 
the natural Hamiltonian of the molecule that provides 
a coupling between the probe and the system qubits of 
the form: Ylj=i Jmc^z 1 ® a< z % ■ This can been seen as 
a controlled rotation of the system qubits about the z 
axis. Therefore, a controlled operation on the system 
qubits can be implemented by allowing the natural evo- 
lution for a time proportional to the desired strength of 
the perturbation. This perturbation can be transformed 
into a rotation about another axis, by sandwiching the 
evolution with two rotations of the target qubits. The 
fidelity decay depends on the eigenvector statistics of the 
system in the eigenbasis of the perturbation F° r 
the large majority of randomly chosen perturbations, a 
regular system will have random statistics in that eigen- 
basis and show an exponential decay at the FGR rate. 
However, for perturbations given by simple functions of a 
system coordinate, only the integrable system will have a 
structured form in the perturbation eigenbasis and hence 
exhibit a fidelity decay with large fluctuations from the 
FGR average. 

Different random maps can be generated by varying 
the pulses and delay times chosen. This allows an averag- 
ing over random maps which is important because of the 
finite size fluctuations in the fidelity decay. The average 
fidelity decay depends only on the relative randomness 
between the system and perturbation. But, the fidelity 
decay for any particular map will fluctuate (an effect in- 
dependent of the initial state which is already averaged 
over by this circuit implementation). These fluctuations 
disappear as the size of the system increases beyond a few 
qubits, but for small systems, this effect is pronounced 
and requires an averaging over different random maps 
to obtain good statistics. To quantify the fluctuations, 
numerical simulations of the fidelity decay experiment 
were performed over 1000 different random maps created 
assuming perfect control. Averages and deviations were 
then calculated to compare with the experimental results. 

The results of Fig. also show how the form of the 
perturbation effects the decay. The pseudo-random map 
shows a universal response, in that the decays are iden- 
tical under different perturbations. On the other hand, 
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FIG. 3: Comparison of fidelity decay curves for regular (a) 
versus chaotic (b) evolution under different forms of pertur- 
bation. The natural evolution of the molecule provided the 
regular system and a pseudo-random map was used to model 
chaotic evolution. The perturbation form was controlled by 
the presence or absence of conjugate rotations on either side 
of the coupling interaction between the probe and system. In 
the case of regular evolution, when the perturbation is a ro- 
tation about the z-axis (A curves), which commutes with a 
system coordinate, the fidelity decay shows substantial fluc- 
tuations away from the average. Changing the perturbation 
to a rotation about the x-axis (■ curves) substantially alters 
the form of the decay. In the chaotic case, the evolution looks 
random in the eigenbasis of either perturbation and the de- 
cays exhibit the universality of the fidelity decay for complex 
dynamics. The FGR result (solid curves) and standard de- 
viations (dashed lines) calculated from numerical simulations 
assuming perfect control are also shown. The inset of figure 
(b) shows the experimental average fidelity decay from twenty 
different random maps. The dotted line show the correspond- 
ing average from the numerical simulations and we can see 
that deviation from the exponential decay is also present in 
the numerical results and well explained by the finite size sat- 
uration level. 



the natural evolution's decay varies wildly. This result 
demonstrates the necessity of choosing the perturbation 
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carefully in the context of distinguishing between regular 
and chaotic evolution: fidelity decay will provide useful 
information only if the appli ed p erturbation commutes 
with the system's coordinate |l3j. 

These experiments also highlight the relevance of these 
techniques to decoherence studies. For this purpose, we 
consider the pseudo-pure qubit to be the system we are 
interested in and the maximally mixed bottom register, 
the environment. The decoherence rate of the system 
is typically governed by some macroscopic parameters 
of the environment, such as its temperature, cutoff fre- 
quency, etc., see e.g. Refs. H3,Q. Recently, the im- 
portance of the dynamics of the environment has been 
expressed by a few authors 0, |2|| , and our experiments 
constitute a direct demonstration of these effects because 
we can directly control the complexity of the environ- 
ment's dynamics. 

Consider a Hamiltonian describing the evolution of the 
probe and environment together as H = He + As ® Be 
where He is the self Hamiltonian of the environment and 
As®Be is the coupling between the environment and the 
probe. We assume an initial product state ps(0) ®ps(0). 
We obtain the state of the probe ps{t) by tracing out 
the environment after a time t. When expressed in the 
cigenbasis of the coupling (A\j) = the diagonal 

elements of ps{t) are equal to those of ps(0), while its 
off-diagonal terms have decayed by a factor 7jfc(i) given 
by (using the Trotter decomposition and S — ► 0) 

—Tr ^(e lS(HE+x i B) y /S (^ e -iS{\j-X k )B e -i5{H B +x i B)y ,S 

which is identical to the trace element of Eq. [21 and so 
was directly measured in our experiments. Hence, the 
decoherence rate in this model depends on how random 
the environment's dynamics appear in the cigenbasis of 
the coupling. Moreover, the coupling strength is modu- 
lated by the eigenvalue difference Xj — Xk , so in the FGR 
regime, the decoherence rate of the (j, k) entry of ps is 
proportional to (Xj — A^) 2 . 

Conclusion — Using our universal control over the qubit 
register we were able to simulate regular and complex 
dynamics for a quantum system and contrast their fi- 
delity decay. Our scheme avoids the scalability limita- 
tions of conventional NMR QIP by relying on neither full 
pseudo-pure state initialization nor full tomography. The 
method is therefore scalable in the absence of noise and 
provides a powerful tool for characterizing the dynamics 
of an unknown system. Moreover, our experiment can 
also be seen as a direct measurement of the decoherence 
rate of a system coupled to an engineered environment; 
wherein, universal control over the environment allows 
for systematic study of the impact of the environment 
dynamics and the coupling to the environment on the 
decoherence rate. As such, it is an experimental demon- 
stration of the importance of the environment dynamics 



in the decoherence rate. 
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